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Fundamental physical phenomena such as laser-induced ionization jTj, driven 
quantum tunneling Faraday waves [HHS], Bogoliubov quasiparticle excitations 
[7J, and the control of new states of matter [SV Hl) rely on time-periodic driving 
of the system. A remarkable property of such driving is that it can indnce the 
localized (bound) states to resonantly couple to the continuum (as illustrated in 
Fig. [^). Therefore experiments that allow for enlightening and controlling the 
mechanisms underlying snch coupling are of paramount importance. We imple¬ 
ment snch an experiment in a special optical fiber characterized by a dispersion 
oscillating along the propagation coordinate, which mimics ’’time” (see Fig. 0>)- 
The quasi-momentum associated with such periodic perturbation is responsible 
for the efficient coupling of energy from the localized wave-packets (solitons 
in anomalous dispersion and shock fronts in normal dispersion) sustained by 
the fiber nonlinearity, into free-running linear dispersive waves (continuum) at 
multiple resonant frequencies. Remarkably, the observed resonances can be ex¬ 
plained by means of a unified approach, regardless of the fact that the localized 
state is a soliton-like pulse or a shock front. 

Our experimental set-up realizes a platform for studying the coupling of localized states of 
a nearly conservative (Hamiltonian) system into radiation modes induced by periodic driv¬ 
ing. The localized states that we exploit are nonlinear (non-spreading) wavepackets of optical 
fibers, namely temporal solitons [131133] and dispersive shock waves (DSWs) [TH - fTU] . They 
share the property to have a well-defined wavenumber and group-velocity which are crucial 
for determining their resonances. Their excitation is readily accessible by operating with 
opposite group-velocity dispersion (GVD, (32 = (Pk/du'^): anomalous {(32 < 0) for solitons in 
order to balance nonlinearity, whereas DSWs emerge from a gradient catastrophe occurring 
in the normal GVD {(32 > 0). Despite such difference, both are essentially described by a 
Hamiltonian model, namely the nonlinear Schrbdinger equation (NLSE) [T31I33] . where the 
propagation distance plays the role of evolution variable (usually time). This allows us to 
introduce periodic driving by employing photonic crystal fibers whose flexibility to engineer 
dispersion is fully exploited to have a longitudinal periodic GVD, thus realizing dispersion 
oscillating fibers (DOFs) [T7].We report clear evidence that this built-in and tailorable pe¬ 
riodicity of the dispersion is responsible for the parametric excitation of radiation modes 
which are amplified out of quantum fluctuations at multiple resonant frequencies (even for a 
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spatially harmonic variation). This thus constitutes a novel implementation of quasi-phase- 
matching [IH]. This is in marked contrast with radiation caused by standard perturbations 
such as third-order (or higher) dispersion [TUl - E^ . which usually feature isolated resonances. 
Our results also clearly show that the parametric excitation of resonances is neither a pre¬ 
rogative of solitons nor of systems with periodic injection or extraction of energy such as 
hber lasers [25] or lumped ampliher links [261 EZ] where dissipative resonances lead to the 
generation of so-called Kelly sidebands. 


Results 


Theory. Perturbation theory allows us to predict the frequency detunings (jJrr of the 
resonant radiations (RR) that can be parametrically excited in a DOF. They turn out to be 
given by the roots of the following expression (see Supplementary information for a detailed 
derivation): 

27r 

D{urr) - kni =-^m, m = 0, ±1, ±2,... (1) 

where £>{ 00 ) = —Akiu + -|- -|- ... is the average linear dispersion in the pump 

reference frame. = d^k/du^ and Aki arises from the deviation of the actual group- 
velocity from the natural one [22]) • In Eq. (18), kni is a well dehned nonlinear wavenumber 


hxed by power. For a bright soliton kni = 7 ^/ 2 , where P and 7 are the soliton peak power 
and the hber nonlinearity, whereas for a shock wave kni = —^£1 being Pb the power of the 


background held over which the RR propagates ra- Equation (jl^ generalizes previously 
proposed formulas [301 El] nnd expresses momentum conservation: it states that the dif¬ 
ference between the momentum of the linear waves and the momentum generated by the 
nonlinear pump must be equal to the virtual momentum carried by the periodic modulation 
of the dispersion. 


We verihed that Eq. (18) describes accurately the parametric excitation of RR in both 


GVD regimes. A clear illustration of this phenomenon is provided in Fig. [^ which shows 
the simulated evolution of a hyperbolic secant pulse (150 fs duration, peak power P) based 
on the NLSE with included periodic GVD and dispersion slope /^s (see methods). Figure 
I3> displays the time domain evolution of a pulse with P = 15 W, corresponding to a nearly 
fundamental soliton launched in a 150 m long DOF of period A = 5 m. Strong radiation 
traveling slower than the soliton becomes visible beyond the activation distance Za ~ 20 m. 
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corresponding to maximum pulse compression. In the spectral domain, the radiation modes 
correspond to several distinct and well dehned frequencies [see Fig. It. c] that agree with 
the prediction based on Eq. (|T^ (vertical green lines). The peak detuned by ~ 9 THz from 
the pump turns out to be the standard RR [m = 0 in Eq. (18)], whereas other twelve peaks 


originate from the periodic perturbation [m 7 ^ 0 in Eq. (18) . 

When pumping in the normal GVD regime, we selected a higher peak power and duration 
(F = 100 W, 280 fs) in order to access wave-breaking, and a shorter modulation period 
(A = 0.5 m) in order to have resonances in a realistic frequency span. In this case, the 
localized state is the shock front that emerges over the pulse leading edge after the breaking 
and activation distance ^ ~ 8 m, which is clearly visible in Fig. mi. Also in this regime, we 
identify a hrst spectral peak at —15 THz in the anomalous GVD region [see Fig. § 3 , 1 ] as 
the standard RR (m = 0 mode) due to /^a, while other Eve peaks are clearly visible. These 
are parametrically excited resonant modes that correspond to m = ±1, ±2, 3 which position 


are also well predicted by Eq. (18)(see vertical green lines). 

We emphasize that in both conhgurations, the m = 0 isolated resonance lies in the oppo¬ 
site GVD regime with respect to the localized wavepackets that generate it, and disappear 
for /^a = 0. Gonversely, the parametric resonances with m 7 ^ 0 lie in both normal and 
anomalous GVD regimes. They also survive the absence of higher-order dispersion (/3n = 0, 
n > 3), and their number is inhnite for a harmonic perturbation. In practice, their actual 
number is limited because their amplitude decreases with the frequency shift and they be¬ 
come weaker than the pulse spectrum envelope and/or the noise background. Note also that 
it is important to take into account the group-velocity deviation Afci [see Fig. [^,d] in the 
D term of Eq. (18) in both conhgurations, in order to accurately predict the actual values 
of the resonances. 

The nature of coupling process is further clarihed by the spectrograms corresponding 
to Fig. which are displayed in Fig. In both cases, the pump wave-packets remain 
clearly localized in both temporal and spectral domains, while shedding energy into radiation 
modes, which disperse in time while remaining at the frequencies predicted by the resonances 


in Eq. (18). We have found that the localized pumps experience temporal and spectral 
breathing with the period of the perturbation and rehll the radiation at each cycle of spectral 
broadening (see multimedia hies). In spite of this transfer of energy the radiation damping 
of the pump is small, thus conhrming the metastable nature of these wave-packets 
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Experiments. We have designed two experiments in order to observe the parametrically 
excited resonances from both solitons and DSWs. The experiment, sketched in Fig. 
and detailed in methods, simply consists of a femtosecond laser pulse launched in DOFs, 
whose parameters are close to those used in the illustration examples in Figs. [2]j^ Due to 
experimental constraints, only the initial pulses were slightly different. In a hrst experiment, 
a 150 fs pulse centered at 1075 nm was launched in a 150 m-long DOF with a modulation 
period of 5 m. Figure shows the longitudinal prohle of the hber diameter measured during 
the draw process (left axis) and the corresponding calculated zero dispersion wavelength 
(ZDW, right axis). At the pump wavelength of 1075 nm, the average dispersion is slightly 
anomalous, so that the pump pulse is able to excite a near-fundamental soliton for peak 
powers in the order of a few tens of watts. An experimental power map, representing the 
output spectrum recorded for increasing pump peak power, is plotted in Fig. It shows 
that, for increasing pump peak power, the spectrum rapidly evolves from a hyperbolic secant 
shape to a much more structured and highly asymmetric one containing more and more sharp 
spectral resonances. More precisely, these discrete spectral sidebands corresponds to the 
parametric excitation of the RR that stems from the periodic variation of the second order 
dispersion. Numerical simulations of the generalized NLSE (see methods and supplementary 
information), without any free parameter, also reproduce these features with an excellent 
quantitative agreement. This can be seen in Fig. B. where the experimental spectrum 
obtained for a pump peak power of 25 W (red line) is compared with the simulated one (blue 
line). In addition, it is worth noting that the theory is quite robust since the observed peaks 


are accurately predicted by solutions of Eq. (18) (green vertical lines) in these experiments. 

The second experiment presented here is devoted to investigating radiating DSWs. The 
DOF is shorter (50 m) as well as the modulation period (0.5 m, see longitudinal prohle 
in Fig. [^). It has been pumped with 280 fs pulses. The pump wavelength was tuned 
to 1037 nm so that it lies in the normal average dispersion region required to generate a 
shock wave from a few hundreds of watts of peak power (see additional experimental details 
in supplementary information). Figure shows the experimental power map. Starting 
at about 50 W, a RR peak is generated across the average zero dispersion wavelength. 
For increasing pump peak powers, additional discrete spectral sidebands corresponding to 
parametrically excited RR appear, similarly to the soliton case. The DOF period A being 10 
times shorter in this case, the spacing between two adjacent peaks is vTO times larger than 
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in the soliton case, as expected from Eq. (18). These results are again in excellent agreement 


with numerical simulations using a generalized NLSE and with Eq. (18), as shown in Fig. 
and in supplementary information. 


Conclusions 

To summarize, we have demonstrated that localized states (either a soliton or a disper¬ 
sive shock wave) can efficiently transfer energy to multiple resonant radiations at different 
frequencies, as a result of the quasi-momentum associated to a DOF. Our experimental re¬ 
sults, supported by numerical simulations and by the perturbative analysis that leads to Eq. 


(18), prove that a DOF is a very simple and highly tailorable platform allowing to study the 


periodicity-induced coupling of nonlinear bound states to the radiation continuum, which is 
a general feature in systems driven by a time-periodic Hamiltonian. The DOF platform can 
also be successfully used to study how this coupling process due to higher order dispersion 
develop in the presence of a train of random solitons arising from spontanous modulation in¬ 
stability or in driven-damped deformations of Hamiltonian systems such as those describing 
passive ring resonators. 

Methods 


Simulation: The results illustrated in Figs. [2]j^ have been obtained from numerical 
integration of the following NLSE for the electric held envelope E = E{z,t) propagating 
along the hber [33] 

dE h(z)d^E , I3,B^E , „ 

Equation is well known to maintain conservative (Hamiltonian) structure despite the 
periodic perturbation embedded in the term (32{z) = (32 + 5(32sm{2'Kz/K) (where 5(32 is the 
perturbation amplitude around the average GVD (32), and the third-order dispersion (3^ ( 
which plays a signihcant role in experiments performed by pumping close to the ZDW, as 
in our case). We have employed the following values of the parameters that arise from hber 
characterization: 7 = 10 (W km)“^, 6(32 = 1-2 ps^/km, (32 = —1.2 ps^/km and (3^, = 0.0716 
ps^/km for the soliton conhguration, whereas (32 = 2.9 ps^/km and (3^ = 0.0645 ps^/km 
for the DSW conhguration. The nonlinear term in Eq. ([^ acts as a self-induced potential 
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that allows for the existence of localized states in the form of bright solitons for /32 < 0 
(anomalous dispersion) or dispersive shock waves for /?2 > 0 (normal GVD), whose leading 
front is, in the unperturbed case, locally a dark soliton [23]. The radiating soliton excited in 
Fig. c is characterized by a soliton number N = slightly higher than 1, so 

that the pulse undergoes temporal compression and hence spectral broadening. It thus acts 
as an effective seed for the parametric excitation of RR. We emphasize that the periodic 
perturbation is nearly resonant with the soliton period (~ 10.5 m), which distinguishes our 
regime from non-radiating guiding-center or average solitons [2H1123], which are normally 
seen when the periodic perturbation is much faster than the soliton period. In the shock 
case, the equivalent quantity N = is much larger than unity {N 3> 1) so that 

the nonlinearity drives the pulse towards the gradient catastrophe that causes the formation 
of the shock front [15] . 

Numerical simulations of Eq. (|^ have been performed by using the split step Fourier 
method with a temporal resolution of 5 fs, 2^^ points and a spatial resolution ranging from 
0.025 m to 0.1 m for the shock and the soliton conhgurations, respectively. 

When directly comparing the numerics with the experimental data (results in Fig. |^, 
we have accounted also for secondary effects in the hber such as losses, Raman scatter¬ 
ing, self-steepening, and fourth-order dispersion by making use of a generalized NLSE (Eq. 
(2.3.36) in Ref. [33], also reported explicitly in the Supplemental information). However, 
we have verihed that the impact of such effects on the pulse evolution over the length and 
temporal scales involved in the experiment is really minor, and most of all does not affect 
the resonances. 

Experiments: Experiments have been performed exploiting a Ti;Sa oscillator delivering 
140 fs full width at half maximum (FWHM) near transform limited pulses. They are sent 
into an optical parametric oscillator (OPO) allowing to generate tunable slightly chirped 
femtosecond pulses. The output beam then passes through a combination of two polarizers 
and two half-wave plates in order to carefully adjust the polarization state and pump power 
simultaneously. It is launched in the DOF with an aspherical lens. The pulses were charac¬ 
terized with a frequency resolved optical gating (FROG) system, before being injected into 
the DOF. For the soliton experiment, the OPO was bypassed so that pulses directly coming 
from the Ti:Sa oscillator were used. They were centred at 1075 nm, and they were measured 
at 150 fs FWHM at the DOF input (i.e. after the launch lens) with a small chirp C = 0.35 
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(as defined in [33]). For the dispersive shock wave experiment, the OPO was used and tuned 
to 1037 nm. The pulses at the DOF input were measured at 280 fs FWHM, with a chirp 
C = 1.24. Spectra out of the DOF were acquired with an optical spectrum analyzer with 
a resolution of 0.2 nm. The output power was measured with a power-meter and the input 
power was deduced knowing the DOF total attenuation. It was cross-checked by cutting the 
DOF at the end of the experiment and measuring the power out of the 0.5 m DOF initial 
section. 
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FIG. 1: Sketch of physically equivalent phenomena, a, coupling of bound states to contin¬ 
uum induced by generic potential which is oscillating in time, b, fiber with oscillating diameter 
(dispersion) which induces wavepackets that are localized in time thanks to a self-induced poten¬ 
tial given by the nonlinearity to couple into the radiation continuum at characteristic resonant 
frequencies. 
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FIG. 2: Parametrically excited multiple resonances at fixed power (numerical model¬ 
ing). Left panels a,b,c refer to the soliton pumped configuration, whereas right panels d,e,f refer 
to the shock configuration. a,d: colormap of the temporal evolution of power (log scale) along 
the fiber; b,e: output spectra superimposed on the dispersion relation D{u}) — kni (red curve). 
The vertical green lines correspond to the resonances determined by the graphical solution of Eq. 


(18), i.e. the crossing between the dispersion curve and the dashed horizontal lines that corre¬ 
sponds to integer multiples of quasi-momentum 27r/A; c,f: colormap evolution of the spectrum 
along the fiber. All the results are obtained by numerical simulation of the NLSE (see methods), 
with periodic GVD and /^s as additional parameter. 
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FIG. 3: Spectrograms corresponding to evolutions in Fig. (numerical modeling), a: 

soliton configuration, b: shock wave configuration. The spectrograms represent the spectra over 
gated time intervals and are computed at the output of each fibers (for the evolution along the 
fibers, see the additional multimedia files). ZDW: zero-dispersion wavelength. 
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FIG. 4: Experimental set-up P: polarizer; HWP: half-wave plate. 
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FIG. 5: Parametrically excited multiple resonances from a fundamental soliton (left 
panels) and from a DSW (right panels): experimental results. a,d: Longitudinal evo¬ 
lution of the DOF diameter (left vertical axis) and consequent zero-dispersion wavelength (ZDW, 
right vertical axis) used for the soliton (a) and dispersive shock waves (d) experiments. b,e: Ex¬ 
perimental power maps showing the development of asymmetric sidebands with increasing pump 
peak power. c,f: Comparison between experimental (red lines) and simulated spectra (blue lines) 
for a pump peak power of 26.4 W in the soliton case (c) and of 234 W in the dispersive shock wave 


case (f). Vertical green lines depict the resonance from Eq. (18). 
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I. SUPPLEMENTARY INFORMATION 


A. Derivation of the resonances 

In order to derive the equation of the resonances [Eq. (1) of the paper], we start from 
the following conservative nonlinear Schrodinger equation (NLSE) 


lOzU + d{idt)u 2 

^-d^U + 7 M = 0 , 

( 2 ) 

diidi) = Y, §(»<)”. 


(3) 


n>2 n>2 


where t stands for time measured in the frame moving at natural group velocity of light 
Vg = dk/d(jj~^ |33], 7 is the nonlinear coefficient, = dP'k/du'^ is the n-th order average 
dispersion, d{idt) and D{uj) are the average dispersion operators in time and frequency 
domain (Fourier transform), respectively. In Eq. (|^ the average group-velocity dispersion 
(GVD) (32 is perturbed by the periodic modulation 6/32{z) and by higher-order dispersive 
effects {(3n, n > 3). 

We consider a localized wave-packet (either a soliton or a shock front) that moves with 
characteristic velocity Vg = 1/Afci, or in other words a solution Us{z,t) = A{T)exp{iqz) 
of the unperturbed Eq. ([^ which has nonlinear wavenumber q and is stationary in the 
reference frame t = t — Akiz. In terms of this new variable Eq. (|^ reads as: 

idzU + d{idr)u - +'j\u\‘^u = 0, ( 4 ) 

d{idr) = ^{idrV — Akiidr, D(uj) = — Akioo. ( 5 ) 

^^ n! ^' n\ 

n>2 n>2 

When the perturbation is effective, the solution can be decomposed into the localized 
wave and radiation modes g{z, r) 

u{z,t) = + g{z,T). ( 6 ) 

Assuming the radiation modes to be weak (\g\ -C |A|), the linearization of Eq. ([^ around 
the localized wave yields 

idzg + d{idr)g - + 7 {u^g* + 2\us\‘^g) = - dH{idr)ug, ( 7 ) 

where dniidr) = J 2 n >3 accounts for the perturbation due to higher-order dispersive 

terms. The RHS of Eq. ([^ represents the driving force for the evolution of the radiation 


16 






modes. Importantly, the effective wavenumber of this forcing is not simply the wavenumber 
of the localized wave-packet but is also affected by the quasi-momentum associated with the 
periodic modulation of the second order dispersion. 

Let us focus first on the free evolution of the system, described by the LHS of Eq. ([^. 
Without loss of generality, the radiation can be searched in the form 

r , ,2 

g{z,t) = exp 


.OJ 


iqz + i— / 6/32{s)ds 

^ ./n 


( 8 ) 


We find the following system that rules the free evolution [i.e. neglecting the forcing term 
corresponding to RHS of Eq. ([7|)] of the dispersive waves 


id. 


+ 


= 0. (9) 


D{oj) — k — q + 2'y\A\‘^ exp 5(32{s)ds\ 

exp [iu"^ S/32(s)ds] - ^D(-u) + k - q + 27 |Ap 

The dispersion relation of the linear waves k = k{u:) can be found by setting the determinant 
of the matrix equal to zero. In terms of the odd and even dispersive contributions {Do{u}) = 
D{u) — D{—uj) and De{uj) = D{uj) + £)(—cn), respectively), we find 


1 

k±(<^) = — ± 2 


£>e(w) - 2^ + 27 |Ap -2q + &i\A\^ 


( 10 ) 


For the soliton, Eq. (10) holds true with A = Q since radiation is not temporally 
overlapped with the soliton (it grows on soliton tails where A exponentially vanishes) and 
q = 7 P /2 [23], being P the peak power of the soliton. Therefore we obtain 


k{uj) = 


( 11 ) 


For the shock wave, the radiation modes are amplified out of noise close to the leading 
edge of the dispersive shock wave. Therefore they propagate over the flat-top background 
Aq (power Ph = l^oP)) which develops as a result of the steepening of the pulse edges [23] . 
which is in turn responsible for the shock formation. In this case we can set q = '^Pb, i.e. 
the wavenumber of the gray soliton associated to the leading edge of the dispersive shock 


wave. In this case Eq. (10) yields 


= £>.(0.) + 47|7l„p 


( 12 ) 


In the shock case the presence of the background generates two symmetric branches of the 
dispersion relation. This fact accounts for the four wave mixing between the dispersive waves 
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and the pump. Usually the amplitude of the two symmetric waves (a and b) are orders of 
magnitude different, so that only one branch of dispersion relation [k^ in our case) turns out 


to be relevant. Under the hypothesis ^ I-Del, we can expand in Eq. (12) the square 
root to obtain, for the relevant branch. 


k{uj) = D{uj) + 'jPb. 


( 13 ) 


The forcing term F arising from the RHS of Eq. ([^ which is effective for the growth of 
the radiation modes with complex amplitudes a and b, turns out to be 


^ S/32{z) 

F = —-— exp 




^ Jo 


d^A — exp 


.00 


i— / S/32{s)ds 


dH{idr)A. (14) 


By considering a periodic d/ 32 {z) with period A, we can expand the perturbation in Fourier 
series of the form 6 fJ 2 {z) = X]/ and consequently expand the exponential in Eq. (14) 

as 

,2 rz 


exp 


.00 


—i— / S[J 2 {s)ds 




( 15 ) 


Therefore we cast Eq. (14) in the form 


F = d^Aj 2 - dH{idr)A ^ d, 

l,m m 


1 


( 16 ) 


which allow to recognize two different driving terms for the growth of radiation modes. The 
hrst one comes from the modulation of the GVD, whereas the second one from higher-order 
dispersive terms. Importantly, even in the case of sinusoidal modulation, the exponential 


term (15) generates an inhnite set of Fourier harmonics. 


Coupling of energy into the radiation modes (dispersive waves) efficiently occurs when 
their wavenumber k{u) equals the wavenumber of the forcing term. This leads to the fol¬ 
lowing resonance condition (quasi-phase-matching) 


27r 

k{u) = m = 0 , ±1, ±2 ,... 

or, equivalently, in the form of Eq. (1) of the paper: 

27r 

D{u) — kni = m = 0, ±1, ±2,. 


(17) 


( 18 ) 


where kni = 7-P/2 and kni = —7-Pb, for the soliton and the shock conhguration, respectively. 
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Equation (jl^ retains its validity for any general dispersion profile. However, in our 
experiments the dispersive operator can be truncated to the hrst correction to GVD, i.e. 
third-order dispersion /ds, whereas all the higher-order dispersive terms can be safely ne¬ 
glected. In this case, Eq. ( [TtI ) can be cast, for the soliton configuration, in the final form 

—UJ + —u - Miu - 7 ^ = ( 19 ) 

whereas for the shock configuration, we obtain 

— Akiuj + 'jPh = -j-m. ( 20 ) 

6 2 A 

B. Properties of the dispersion oscillating fibers 

Two different DOFs with different lengths and modulation periods have been fabricated 
for the experiments. The fiber used for the soliton experiment, labelled DOF^l here, is 
150 m long and has a modulation period A = 5 m. Its outer diameter, which longitudinal 
evolution is displayed in Fig. 5a of the manuscript, oscillates between 110 and 123 pm. 
The inset in Fig. shows a scanning electron microscope (SEM) image of the DOF cross 
section. It has two bigger holes around the core allowing to increase its form birefringence and 
to ensure a polarization-maintaining behaviour. Figure shows the full dispersion curve 
simulated with a commercial finite-element mode solver for the maximum and minimum 
diameters (red and blue lines respectively), for the neutral axis excited in experiments. The 
black line represents the average dispersion over the whole DOF length. The average zero- 
dispersion wavelength is located at 1064.5 nm. Its attenuation was measured to be 0.75 
dB at 1064 nm and its nonlinear parameter was calculated to be 7 = 10 (W km)“^ at this 
wavelength. The DOF used for the dispersive shock wave experiment, labelled DOF^2, 
is 50 m long and has a modulation period A = 0.5 m. Its outer diameter (Fig. 5d of the 
manuscript) oscillates between 110 and 122 pm. It is based on the same design than DOF^^/:! 
although the geometrical parameters of the cross-section are slightly different. Figure 
shows the simulated dispersion curves for the maximum and minimum diameters (red and 
blue lines respectively) as well as the average one (black line). The average zero-dispersion 
wavelength is located at 1062.5 nm. The overall attenuation is 0.35 dB at 1064 nm, and its 
nonlinear parameter is 7 = 10 (W km)“^ at this wavelength. 
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Parameter 


Soliton 


Shock 


/32 (ps^/km) 

-1.2 

2.9 

5(^2 (ps^/km) 

1.2 

1.2 

ds (ps^/km) 

0.0716 

0.0645 

(3 4 (ps^/km) 

-1.1 • 10“^ 

-1.1 • 10“^ 

7 (W“^km“^) 

10 

10 

a (dB/km) 

5 

7 

A (m) 

5 

0.5 

L (m) 

150 

50 


TABLE I: Parameters of the fibers employed in the numerical simulations 

C. Numerical simulations and fiber parameters 


The parameters of the fibers that arise from fiber characterization and have been used 
in numerical simulations are listed in table 1. In the illustrating example reported in Figs. 
2 and 3 of the manuscript, only the elements which are essential for the description of the 
basic phenomenon have been considered. In this case the simulations has been performed 
by integrating the nonlinear Schrodinger equation (NLSE) reported in the method section, 
i.e. the basic NLSE with additional perturbations due to third-order dispersion and the 
periodic GVD 6 / 32 {z). For simplicity, the initial conditions have the following hyperbolic 
secant shapes E{t,z = 0) = A/^sech with Tq = 150 fs and Pq = 15 W, for the 

soliton configuration, and Tq = 280 and Pq = 100 W for the shock wave configuration. The 
spectrograms [Fig. 3 in the paper] have been calculated by using a Gaussian pulse as a gate 
with duration of 1.6 ps. 

When comparing directly with the experimental results [see Fig. 4c,f in the paper] we 
make use of the following extended NLSE [3^, which accounts also additional effects such 
as higher order dispersion effects, Raman effect, self-steepening, and fiber losses 

dE{z,t) _ .( 32 {z) d^E{z,t) / 33 d^E{z,t) 64 d^E{z,t) «p. 

dz ~ ^ 2 dP ^ 6 dP ^^24 dP 2 ^ 

E{z,t) j R{t')\E{zp-t')\dt' 

with Pn the nth order dispersion terms, a the linear losses, Ts = l/cjp with ujp the central 
pulsation of the pulse, R{t) the full nonlinear response function that includes the instanta- 





1 • ^ 


X 
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neous (Kerr) and delayed (Raman [M]) contributions with fractional weights fKerr = 0.82 
and fRaman = 0.18, respectively [33]. Here P 2 {z) = P 2 + A/32sin(27r2;/A) and we checked 
through numerical simulations that all other parameters can be assumed to be constant 
along the hber length. Indeed, the modulation of the nonlinear coefficient for instance is 
about 10 %, which is about one order of magnitude lower than the /32 one [35] . 

In this case we also use initial conditions that accurately describe the pulses injected in the 
hber, which has been experimentally characterized by means of the frequency resolved opti¬ 
cal gating (FROG) system. In particular the best ht with FROG data gives slightly chirped 


input pulses of the following form: E{t, z = 0) = A/^sech 


1.76t 

To 


exp( 


-iCt^ ' 


2T§ 


for the soliton 


conhguration with Tq = 150 fs and C = 0.35, and E{t,z = 0) = \/^6xp 


-1.665(l+iC)t^ 


2T2 


for 


the shock conhguration with Tq = 280 fs and C = 1.24. The outcome of the numerical inte¬ 


gration of Eq. (21) with such initial conditions are directly compared with the experimental 
results in Fig. 4c,f of the paper. 

Using the realistic simulation parameters given above, we have also simulated the spectral 
output against pump peak power in both the soliton and dispersive shock wave conhgura¬ 
tion, respectively. The results, displayed in Fig. 01 ,b are directly comparable with the 
experimental maps reported in Figs. 5b,e of the paper. In both cases the quantitative 
agreement with experimental data is excellent and conhrms the progressive excitation of 
multiple resonances induced the periodic driving with increasing peak power. 

Finally, we emphasize that we have compared the experimental results with the simu¬ 


lations of the full model [Eq. (21)] in order to have a better accuracy over all the details. 
However, we have verihed that no signihcant diherence arises when integrating the simpler 
NLSE reported in methods, since the additional terms (steepening, Raman ehect, losses, 
fourth-order dispersion) are indeed negligible in our regime. In particular the parametric 
excitation of the resonances that we observe in the experiments and in the numerics are 
indeed quantitatively and accurately explained on the basis of this simpler NLSE. 
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FIG. 6: Simulated dispersion curves of the DOFs used for the soliton (a, DOF^l) and disper¬ 
sive shock wave (b, DOF^2) experiments. Red and blue lines correspond to the maximum and 
minimum diameters, respectively. Black lines correspond to the average dispersion curves over 
the whole DOF lengths. Insets show SEM images of the DOF cross sections at the maximum 
diameters. 
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FIG. 7: Colormap of the output spectrum against the input power obtained by means of numerical 


integration of Eq. (21): a. soliton configuration; b. shock confignration. The dashed vertical lines 
indicates the average zero dispersion wavelength (ZDW). The maps in panels a and b are directly 
comparable with the experimental results in Figs. 5b and 5e of the paper, respectively. 
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